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Abstract: 

The formation and dissociation of bosonic molecules in an optical lattice, formed by spin-1/2 
fermionic atoms, is considered in the presence of an attractive nearest-neighbor interaction. A 
mean- field approximation reveals three different phases at zero temperature: an empty, a con- 
densate and a Mott-insulating phase. The density of fermionic atoms and the density of bosonic 
molecules indicate a characteristic behavior with respect to the interaction strength that distin- 
guishes between a dilute and a dense regime. In particular, the attractive interaction favors the 
formation of molecules in the dilute regime and the dissociation of atoms in the dense regime. 

1 Introduction 

Feshbach resonances provide a powerful tool to vary the interaction strength of atomic gases 
Physical effects, like the formation of molecules, Bose-Einstein condensation or Cooper pairing in 
fermionic systems depend strongly on the strength of the interparticle interaction |21 El |^ E] . 
A model for spin-1/2 fermions with attractive nearest-neighbor interaction is studied to describe 
the formation and dissociation of bosonic molecules in a grand-canonical gas of fermionic atoms. 
The model depends on three tunable parameters, the fermionic tunneling rate t, the fugacity C 
of the grand-canonical gas, and the molecular tunneling rate J. J also controls the attractive 
nearest-neighbor interaction of the model. 

A difference in comparison with previous works [3 El El ^3 is that we study the T = 
properties of a Fermi gas in an optical lattice. A consequence of the latter is that a Mott-insulating 
state can be identified, whereas T = avoids thermal fluctuations of the system. In this case a 
mean- field approximation should give reasonable results, provided that quantum fluctuation are 
not too strong. 

A qualitative picture of the T ~ phase diagram can be given by an estimate of the energies. 
There is a single-particle potential due to the chemical potential of the fermionic atoms. It is 
assumed that the chemical potential of the molecules is the same as that of the constituting 
fermionic atoms, i.e. the effective chemical potential of a molecule is twice as that of a single 
fermionic atom. If this potential is sufficiently negative the optical lattice is empty, since particles 
cannot overcome the potential barrier. The kinetic energy, however, represented by the tunneling 
rates t and J, competes with the chemical potential and favors the creation of a Fermi gas and/or a 
molecular Bose gas. The interaction plays a crucial role in this regime: In addition to the formation 
of (local) bosonic molecules, the attractive interaction between the fermionic atoms may lead to 
the formation of Cooper pairs, which can be considered as "non-local molecules" . For sufficiently 
large chemical potential the optical lattice will be filled with two atoms per lattice site, which is a 
Mott-insulating state. This Mott state differs from the Mott insulator of a half-filled Fermi lattice 
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gas with repulsive interaction (Hubbard model ^3]), where the interaction creates a ground state 
with one fermion per lattice site. In other words, the Mott insulator in the system with attractive 
interaction is a special kind of insulator that is created by the Pauli exclusion, and not by a finite 
repulsive interaction, in contrast to the Hubbard model at half filling. 

A mean-field theory is applied to discuss a Mott insulator and a condensed phase. Depending on 
the tunable parameters of the model, the density of dissociated atoms and the density of molecules 
in the condensed phase are investigated. It will be shown that the density of dissociated atoms 
increases with the total number of particles, reaches a maximum and decreases, whereas the density 
of molecules is a monotoneously increasing function. 

2 Model: Hamiltonian and Functional Integral 

A gas of spin-1/2 fermions in an optical lattice is considered. Using fermionic creation (annihilation) 
operators cj^. (cra) for spin <t and at site r (a minimum in the optical lattice), the Hamiltonian of 
the Fermi gas is 

<r,r'> (T=t,X <r,r'> r (T=T.i 

<r,r' > refers to nearest-neighbor sites r and r'. A chemical potential term for molecules 

Et t 

r 

has been neglected because it leads only to a shift of J in terms of the subsequent mean-field 
calculation. Similar Hamiltonians were considered in a number of papers [Sjl^^]. The first 
term describes tunneling of individual fermions in the optical lattice with rate t, the second term 
tunneling of local fermion pairs (i.e. bosonic molecules) between nearest-neighbor sites with rate 
J. It should be noticed that the latter is also responsible for an attractive interaction between 
fermions with different spin (J > 0). However, there is no local (diagonal) interaction, all the 
fermionic interaction is carried by the nearest-neighbor term. The bosonic molecules experience 
a local repulsive (hard-core) interaction, since their constituent atoms are fermions and obey the 
Pauli exclusion. For the limiting case i — only bosonic molecules can appear, for J = only 
non-interacting fermionic atoms. The chemical potential fi controls the number of particles in a 
grand-canonical ensemble. The latter is given by the partition function 

Z = Tre-f'''. 

Space-time correlations of fermions are described by the Green's function 



Gr,t,<7;r' ,0,a' — '^'^^ 



The partition function can also be written in terms of a Grassmann integral [131 115| as 

Z = y e-^V[^], 

where the action S for spin-1/2 fermions with attractive interaction, related to the Hamiltonian in 
Eq. iP), reads 



E 



r,r' ,t 



J ^2 ^r,t,l'0r',i+l,l^r,t,2'0r',t+l,2 (2) 



r.r't 
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with space-time coordinates x = {r,t). Wr.r' is l/2d for nearest-neighbor sites on a d-dimensional 
cubic optical lattice and zero otherwise, and ~ \s, the fugacity. The Green's function then 
reads 



After renaming ip hy a time shift i — > t — 1 

A,t,] A,t-i.j {j = 1,2) (3) 

the action reads 

S' = ^{ipx.ldftijx.l + 1px.29'[i>x,2) - ^ iC^r.r' + tWr,r'){lpr.t,l1pr' ,t,l + A.t,2A',t,2) 



-J ^ 1pr,t,l1pr' ,t,li>r,t,2lpr' ,t,2- 
r,r' ,t 



The quartic interaction term is now diagonal with respect to time. It can be decoupled by two 
complex Hubbard- Stratonovich fields (j)x, Xx CSl- The linear combination i(l)+x couples to fermions 
like 

i'ifpx + Xx)ipx,i4'x.2 + h.c. (4) 

This coupling is similar to the atom- molecule coupling discussed in Refs. ^2 ^] ■ There is a 
difference, however, by the fact that we started from the fermionic model and derived an effective 
electron-boson model, whereas the other authors started from an electron-boson model and derived 
the effective fermion model. In both cases an effective coupling constants can be evaluated from 
the original model, an effective fermion-fermion coupling in Ref. and an effective fermion- 
boson coupling in our model. Although this is an interesting problem it will not be pursued in this 
paper. Instead, the fermions are integrated out in since they appear only in a quadratic form in 
the action. This step provides an effective model only for bosons: the integration gives a fermion 
determinant and the resulting fuctional Z depends only on the complex fields (f) and x with the 
effective action 

5e// = (0, v-'^) -t- x) - i log(det A) (5) 

with the antisymmetric space-time matrix 

i<f> + X C + tw-df 

-i(t)-x C + tw-dJ' 

-Q-tw + dt _0 i(l) + X 

-C-tw + dt -i4>~X 



A 



with = {w + 21)-'^ /J. 

(j) corresponds to the conventional BCS field in the case of a local interaction. The additional 
field X is necessary for the nearest-neighbor interaction to ensure that the Hubbard- Stratonovich 
decoupling is well-defined |15) . It will be discussed subsequently that the mean- field approximation 
yields a linear relationship between these two fields. 

2.1 Atomic and Molecular Densities 

The densities of atoms and molecules can be measured as expectation values of the Grassmann 
fields. For this purpose the original fields are used, i.e., the fields before the time shift in Eq. (j^J, 
to write 

.e-'^V[il)\/Z. 
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The action S of Eq. can be expanded in Z around 



X 

to obtain 



l>0 

This expansion can be viewed as an expansion in terms of world hnes in a space-time lattice. There 
are two types of world lines, individual fermion lines with spin a and molecular world lines, as 
shown in Fig. 1. For a given point in space and time x the contribution from the polynomial 
(So — Sy is either a factor 1 (i.e. no contribution from {Sq — S*)' at this point), a factor 4'x,ai'x,a: 
or a factor ipx,iipx,itpx,2''px,2- To measure the probability of the appearence of these factors in Z, 
an expectation value with respect to the Grassmann field can be introduced. In particular, the 
density of dissociated atoms reads 

nf,x = -{ll^x,ltpx,l + 1px,2ll>x,2 + '2i>x,lll>x,ltpx,2i>x,2) 

and the density of molecules 

(1 + tpx.l^x.l + ^x.2i'x,2 + 1px,l'4^x.l1px,2i}x,2)- 

The truncated expectation value 

Cl2 = {lpx,l1px,li>x,2tpx,2) - {^x,li>x,l){tpx,2i>x,2) 

gives 

"■/,x = -(1 + 'Ipx,l1px,l){i>x,2-il^x,2) - (1 + ^x,2lpx,2){lpx,l'lpx,l) " 2Ci2 (6) 

and 

n-m.x = (1 + i>x,li^x.l){l + 1px.2'4'x,2) + C'12. (7) 

It will be seen below that C12 vanishes in mean-field approximation. 

In the Mott-insulating phase a fully occupied lattice is expected, i.e., two fermionic atoms (= 
one molecule) per site, as the only commensurate state, since there is no repulsive interaction which 
could maintain a commensurate state. Any other groundstate of the system is a condensate, except 
for the empty lattice. The remaining question is how the condensate varies with the total number 
of atoms and molecules in the system. It is obvious that the total density of particles increases 
monotoneously with the fugacity. However, it is less clear how the formation of molecules or the 
dissociation of atoms is affected by an increasing fugacity. Moreover, a naive picture suggests that 
the density of molecules increases monotoneously with an increasing attractive coupling of the 
fermionic atoms J. It will be shown in terms of a mean-field approximation that this is not the 
case. 

3 The Saddle-point Equation 

The saddle-point condition for the effective action in Eq. (|3J) is 

5S,ff = 0. (8) 

This gives a nonlinear difference equation. A simple ansatz for its solution is given by uniform 
fields 4> and x (mean- field approximation). Then Eq. ^ reads 

= 3JG(0-ix) 
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with the integral 

G = (10) 

J~i V(-«X + '/')^(-iX + (t>y + 2(s^ + l)(-ix + 0)(-iX + 0) + (s^ - 1)^ 

and s = + tx. p{x) is the density of states for the nearest-neighbor tunneling term w. 
There is a trivial solution (f> = x = and a nontrivial solution 

G-l/J, X- -2*0/3. (11) 

This is the same mean-field equation as in the BCS theory, if J is considered as the coupling 
constant of the fermions. 

The solution in Eq. 111|) implies a non-negative value for 

i-tX + <t>)i~^X + ^)^^^/9=\(t>\V^- 

Then G in Eq. ^ reads 

r dx (12) 

i-1 v/|0lV81 + 2(s2 + l)|0|2/9 + (s2 _ 1)2 

which descreases monotoneously with |0|2. If (s2 _ i)2 becomes inside the interval of integration 
— 1 < cc < 1 the integral diverges for \(j)\ 0. Therefore, there is a non-zero solution of (j) for any 
value of J > 0. On the other hand, if {s"^ — 1)2 does not become inside the interval of integration, 
there is a non-zero solution only for sufficiently large values of J (cf. Fig. 2). 

The densities can be evaluated within the saddle-point approximation. A straightforward cal- 
culation gives for the expectation values, used in nf_x and Um^x, the expression 

{^x dx ,) = -- + - t |0|V9 + s'-l 

^ ''^ 2 2 v/|<^IV81 + 2(s2 -f l)|0|2/9 + (s2 - 1)2^^ ^ ^ ^ 

and 



4 Results 

Using the mean-field approximation of the previous section the following quantities are evaluated: 
the order parameter |(/)|2 of the condensate, the T = phase diagram, and n/(|(^|2,^), nm(|0P,C)- 

4.1 0=0 

In the case of a vanishing order parameter (i.e. (/> = 0) only independent fermions are described 
by the mean-field approach. Then the diagonal element in Eq. (|13|l is 

{-^xAx,]) = - J 6(1 - (C + txf)p{x)dx. 

If s2 < 1 for the entire interval of integration, i.e., 

-l+i<C< 1-^ 

both densities vanish: rif = rim = 0. On the other hand, if s2 > i for the entire interval of 
integration, i.e., 

C > 1 + t 

one obtains = 0, rim = 1- Thus all lattice sites are occupied by bosonic molecules. There is 
also an intermediate regime for 1 — t<(^<l+t where n^, rim 0. Details of this behavior are 
shown in Figs. 2 and 3. 
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4.2 0^0 

For the beginning the special case t — (molecules cannot dissociate) is considered. Then the 
integrand in Eq. H13|l is constant with respect to the integration variable x and yields 

^ -1/2+ i- WV9 + C^-1 



2VI*IV81 + 2K» + l)l0P/9+K»-l)' 



The saddle-point solution 
exists for 



|</>P = 9( + J2 - 1 - C^) 

f V]~J < C < \/ TTj if < J < 1 
1 < C < \/TTJ if J > 1 



The corresponding phase diagram is shown in Fig. 4. This result yields 



= -1/2 - 1/ J + + 1/4. 

In the general case with dissociated atoms (i.e. for t > 0), the density of states in G can be 
approximated by a constant as p{x) — 1/2. Moreover, for small J the integral in Eq. (|12|l can be 
easily performed, giving the saddle-point solution of the order parameter 

~ 4Pe-4*/-^ il-t<C<l + t). 

The model parameters C, J, and t control the densities. However, in a realistic situation the 
order parameter and the densities n/, Um will be measured. Therefore, the densities for a given 
value of\(j)\ are plotted as functions of the fugacity (, where increasing of |i?!>| means an intereasing 
interaction parameter J. The result is shown in Figs. 2 and 3: (1) The density of dissociated 
atoms Uf has a maximum. This maximum is at C = 1 for cj) — and moves to lower values of 
C with increasing order parameter values. The maximum value oi nj itself is 0.5 but decreases 
for |(/)| > 1. (2) The density of molecules rim always increases with But as a function of \ip\ it 
increases (decreases) below (above) C « 1. This result indicates that strong interaction favors the 
formation of molecules in a dilute system and favors the dissociation of atoms in a dense system. 
The latter can be understood as a special kind of frustration effect because there many ways for 
a pair of fermionic atoms to form a molecule. (3) The total density rif -\- rim increases always 
monotoneously with (. 



5 Conclusions 

An attractively interacting Fermi gas in an optical lattice was treated in mean-field approximation. 
At zero temperature three different phases were found: an empty phase, a Mott-insulating phase 
and a condensed phase. The latter is characterized by a non-vanishing order parameter. The 
fermions appear in this phase as a mixture of local pairs (molecules) and extended pairs (dissociated 
atoms) . The densities of these two types of fermionic pairs have a characteristic behavior for the 
crossover from weak to strong attractive interaction. This is indicated by (1) a maximum of the 
density of dissociated fermionic atoms and (2) support for the formation of molecules (dissociation 
of atoms) in the dilute (dense) system by the attractive interaction. 
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Figure 1: Statistics of world lines in a pure system of bosonic molecules (a) and in a mixture of 
fermionic atoms and bosonic molecules (b). Full and dashed lines distinguish the fermionic spin "f 
and [. 
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Figure 2: The density of individual fermions (dissociated atoms) (first figure) and molecules (second 
figure) for different order parameter values (j> = 0, = 0.1, = 0.5, = 1, and = 2, all 
with atomic tunneling rate t = 0.25. 
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Figure 3: The densities as a function of the order parameter. Numbers refer to the value of the 
fugacity ( and letters to fermionic atoms (f) and molecules (m). 




Figure 4: T = phase diagram for atomic tunneling rate t = 0: fugacity C vs. molecular tunneling 
rate J. The indicated densities n are molecular densities Um- The condensed phase with |^| > 
becomes wider when t > 0. In particular, the boundaries of this phase at J = are 1 — t and 1 + t. 
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